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Abstract 

In higher-dimensional supersymmetric theories gauge coupUngs of the effective 
four-dimensional theory are determined by expectation values of scalar fields. We 
find that at temperatures above a critical temperature T^, which depends on the 
supersymmetry breaking mass scales, gauge couplings decrease like T~", a > 1. 
This has important cosmological consequences. In particular it leads to a relic 
gravitino density which becomes independent of the reheating temperature for 
Tr > T*. For small gravitino masses, <^ rrig, the mass density of stable 

gravitinos is essentially determined by the gluino mass. The observed value of cold 
dark matter, Qcdm^^ ~ 0.1, is obtained for gluino masses nig = O(lTeV). 



In higher-dimensional supersymmetric theories [1], where the standard model 
emerges as low-energy effective theory, gauge and Yukawa couplings are determined by 
expectation values of gauge singlet 'moduli' fields. In a cosmological context, this implies 
that generically all couplings depend on the parameters of the cosmological evolution, 
such as the Hubble parameter, temperature, or the cosmological constant. 

In the following we study the dependence of gauge couplings on temperature. As we 
shall see, this has important consequences for the production of gravitinos in the early 
universe. 'Vacuum alignment' at high temperatures causes a power-like decrease of gauge 
couplings. This then leads to a relic gravitino density which becomes independent of the 
reheating temperature Tr above a critical temperature T*. 

As a specific example, consider gaugino mediation [2, 3] which is an attractive mech- 
anism to generate a realistic mass spectrum of gauginos, higgsinos and scalar quarks and 
leptons in the supersymmetric standard model. The source of supersymmetry breaking 
is the vacuum expectation value of a gauge singlet chiral superfield 

(S) ^So + 99Fs , (1) 

which is localized on a four-dimensional (4d) brane embedded in D-dimensional space 
time. The coupling to bulk gauge fields, expressed in terms oi Ad N — 1 superfields, is 
given by 

Id ^ I d^x d^-^y d^O | -^WW 

+ 5i'^-^){y-ys)^SW'^W'^+..^ + h.C., (2) 

where W"' is the supersymmetric field strength. M is a mass scale in the range between 
the compactification scale and the D-dimensional Planck mass, 

^ < M < Md < My> . (3) 
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Here V = J d^ *y is the volume of the compact dimensions, M^, = (VMj^ ^) ^^"^Mp 
and Mp = l/(87rGjv)^/^ = 2.4 x 10^^ GcV is the 4d Planck mass. For instance, with 
l/1/V(o-4) ^ Mgvt = 2 X 10^6 GeV one obtains Md ^ 2 x 10^^ GeV in the case D = 6. 

Inserting the expectation value (1) in the action (2) one obtains for the 4d gauge 
coupling and for the gaugino mass. 



^ + t = A, (4) 



= y^, (5) 
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where 0o = Re-S'o. For the SU(3) gauge couphng of the standard model one has Qodji) > 
g'Q(MGUT) — |- The gravitino mass is given by 

m3/2 = r,^ , (6) 

where 1] > l/VS. The smallest gravitino mass is obtained if Fs is the only source of 
supersymmetry breaking, which is the case in gaugino mediation. The gravitino mass is 
then always smaller than the gaugino mass rrig, 



> f ^ = ^g'oiVME-r^my, > , (7) 

since the volume enhancement factor p = {V M^^'^)^/'^ is larger than 2/{\/3gQ) < 4/\/3. 
For instance, in £) = 6 one has p ~ 10. 

The 4d effective action for the zero modes contains a coupling of the scalar field 
to the supersymmetric gauge kinetic term, 

h ^ j d^x |-^F;^F«'^'' - iAV'^(D^A)" - ]^m-g (A'^A" + A" A") 

+ ("i^M.^""'' - i^"^"(^MA)") +■■■}; (8) 

here F° is the field strength of the vector potential A", and A" denotes the gaugino. At 
finite temperature the gauge kinetic term acquires an expectation value which leads to 
a force on the scalar field (f). This expectation value can be easily calculated by making 
use of the anomalous divergence of the supercurrent [4] , 

5"Jad = \^D^W''W^ , (9) 
which contains the trace anomaly of the energy momentum tensor, 

T^'^^-2^^^(-JF^''F''^''' -iX'a^^D^xY^ , (10) 



90 

where /3(g'o) is the usual /3-function of the gauge coupling. 

The thermal average of the energy momentum tensor is determined by energy density 
and pressure, 

{T^,)t = £ - 3P , (11) 

which are related by 

e^-P + Ts = -P + T^. (12) 
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The pressure has been calculated in perturbation theory for a gauge theory with fermions 
in the fundamental representation [5]. Correcting for the colour charge of the gauginos 
one obtains for a pure supersymmetric gauge theory, 

P = (ao - a^gliT) + ..)t\ (13) 

with 

TT^ 1 
«o = i^riA , 02 = -^AUA ■ (14) 

Here T4 is the Dynkin index of the adjoint representation and Ua = dim G, i.e. the 
number of gluons. For SU(iV) one has Ta = N and ua = N"^ — 1- From eqs. (10) - (13) 
one obtains for the thermal expectation value of the gauge kinetic term, 



__^a^^aM. - iX\a>^D,\r^^ = a^g'^T' . (15) 

Note that the sign of the expectation value is positive and that there is no dependence 
on the /3-f unction. Because of the anomaly one no longer has P — e/3. 

The mass of a chiral superfield, whose vacuum expectation value breaks supersymme- 
try, is generally controlled by the supersymmetry breaking mass scale, i.e., oc 7713/2. 
Small fluctuations around the minimum are then described by the lagrangian (cf. (8), 
(15)), 

jC^l(d<l>r-^-ml/,<l>' + a,gtT'^. (16) 

Hence, the thermal fluctuations of gauge bosons and gauginos induce a negative linear 
term in the effective potential for (f). In many models the parameter ^ is 0{1). 

The negative linear term in the effective potential leads to an increase of the field 0. 
Its equilibrium value at finite temperature is given by 

, «29o" T' 

Note that the fluctuations of are not in thermal equilibrium and that does not acquire 
a thermal mass. According to (4) the shift in (f) changes the gauge coupling to g{(j)T), 

1 + ^ = ^— (IF,) 
gl M g\<i>T) ■ ^ ' 

This change of the gauge coupling becomes significant at a temperature where (pr/M 
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Here we have assumed that Fs does not depend on temperature, as in the Polonyi 
model. Using eqs. (5) and (6) the mass scale M can be expressed in terms of gaugino 
and gravitino masses, which yields 

1/4 /^2, V2 



— ■ (20) 

Extrapolating eqs. (17) and (18) to temperatures larger than leads to a rapid decrease 
of the gauge couphng as 5'^(0r) oc 1/T^. 

However, at large values of (f)T/M the effective lagrangian (16) is no longer appropri- 
ate. First, the decrease of the gauge coupling reduces the force of the thermal bath on 
the field (j). This backreaction can be taken into account by using as effective potential 
the free energy density of the thermal system evaluated with the field-dependent gauge 
coupling, 

/ = -P = (-ao + a29\T, </>) + ...), (21) 

where g{T, 0) has to be determined from the equations of motion. Second, for large values 
of 0, higher powers of (f)/M have to be taken into account. This leads to the effective 
lagrangian 

^ = J (50)' - ^^3/2^(0) - «2/(T, 0)r^ , (22) 



l..,.o 1 

2 

where 



k{(f)) replaces the linear term 0/M in eq. (8). The equilibrium value of is now determined 
by the equation 



h/2 



For small values of one recovers cq. (17). Neglecting corrections 0{(f)/M) for h{(f)) and 
/c(0), keeping only the effect of the back reaction, one obtains at large temperatures 0t oc 
T^/^ and correspondingly for the gauge coupling g'^(T,(f)T) oc T~^/^. This decrease with 
temperature is much weaker than the T~"^ fall-off obtained in the linear approximation. 
We expect that the true decrease, which is determined by the back reaction together 
with the behaviour of h and k at large values of 0, lies somewhere in between. 

The time evolution of the field is determined by the equation of motion 

4> + 3H<j)+ -mloh'U) o A;'(0)r^ = , (26) 
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where H is the Hubble parameter. For H > 7713/2 the motion is damped whereas for 
H < 7713/2 the field (f) oscillates. During the period of reheating the Hubble parameter 
generally depends not only on the thermal bath, but also on the time evolution of other 
fields, in particular the inflaton. The detailed analysis of the time evolution of (j) is 
beyond the scope of this paper. In the following we shall assume that at the end of 
reheating thermal equilibrium is achieved and that, to good approximation, is close to 
its equilibrium value (pT- 

The power-like fall-off of gauge couplings at high temperature, oc T"" with a > 1, 
has important cosmological implications. An immediate consequence is that one loses 
thermal equilibrium at a temperature Teq much below the unification scale Mgut- For 
instance, for a = 2, rrig ~ 1 TeV and 7773/2 ~ 100 GeV, one obtains r(Teq) ~ H{T^ at 
Tcq ~ {rn^i^M'^/mgY/'^ ~ 10^^ GeV. The decrease of the gauge coupling also crucially 
affects the production of gravitinos after infiation [6] which we now discuss. 

The thermal production of gravitinos by gluons, gluinos, quarks and squarks is gov- 
erned by the Boltzmann equations. The collision term has been calculated to leading 
order in the gauge coupling. For the gauge group SU(A'"), with 277/ chiral multiplets in 
the fundamental representation, one has [7] , 

'^''^'^ +mn^,2 = C3/2{T,<P) 



dt 



where 



T{T) = (in (^^——] + 0.3224^ {N + 77/) + 0.578177/ , (28) 



with the thermal gluon mass 



gluon 



ml^o.{T)^^j{N + nf)T'. (29) 



For the gauge coupling we use g{T,(pT)-i except in case of the gluon mass which enters 
only logarithmically. 

In the supersymmetric standard model gravitino production is dominated by QCD, 
the strong interactions, where we have = 3 and 77/ = 6. If the gravitino is the LSP and 
the GUT relations for gaugino masses hold, one has 77^3/2 <C rrig. Integrating eq. (27) up 
to a reheating temperature Tji> T^, assuming a power decrease of the gauge coupling. 
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one obtains a number density to entropy density ratio of gravitinos which is independent 



nz/2 



s To s(T,)H(T,y 

Here Tq is the present temperature, s = (27r^/45)5'* (T)T^ is the entropy density, with 
g*{T^) = 915/4 in the supersymmetric standard model, and 

for q; = 1 . . .4. Inserting the expression for the coUision term in eq. (27) one finds for the 
energy density to entropy density ratio of gravitinos (Tr > T^), 

Hn 135v/l0C(3) iv^-ir.m?(r.) ,,,, 

At temperatures Tr much larger than also contributions involving Yukawa interactions 
may become important, which remains to be studied. 

One can now insert the relation (20) between the temperature T* and gluino and 
gravitino masses into eq. (33), which yields the result (Tr > T*), 



P3/2 



135V5C(3) - 1 Uf(p)\ I e. ^ 

Here we have used the gluino mass at a scale /i as parameter, and !p{T^) = 
T{T^)gy^ {T^) / gl{ix) is a factor 0{1) which takes gauge couplings and their running 
into account. Remarkably, in P3/2/S the dependence on the gravitino mass has dropped 
out. For the dominant QCD contribution N — 3 and 02 = 3/8 (cf. (14)). Dividing by the 
critical density Pait/s — 3.65/i^ x 10~^ GeV [8] one finally obtains {Tr > T*), 



^>V.*^ = 0-l-f^^^)"'(|V^(./m. (35) 



For gaugino mediation one has ^/rj'^ — 0{1); in the temperature range = 
10^ • • • 10^^ GeV we estimate I(^a)^{T*) = 0.5 ... 2. It is then very astonishing how close 
the obtained value for ^13/2/1^ is to the observed one for cold dark matter for gluino 
masses 0{1 TeV). The WMAP collaboration recently obtained {2a error), Qcdm^^ = 
(Qm — ^b)h'^ — 0.113l'lo;oi8 [9]- The relic gravitino density $13/2/1^ is shown in fig. 1 as 
function of the reheating temperature Tr for different values of rrig and 771.3/2. At Tr ~ 
the density reaches a plateau whose value is essentially independent of Tr and 777.3/2. The 
figure clearly shows the scaling oc 7773/2/^/777^. 
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Figure 1: Relic gravitino density ^3/2/1^ as function of the reheating temperature Tr 
for different gravitino and gluino masses: 1713/2 = 20 GeV with rrig = 1.5 TeV (dashed 
hue), m-g = 1.0 TcV (full hue), m-g = 0.5 TeV (dotted line), and ms/2 = 200 MeV 
with rrig = 1.0 TeV (dashed-dotted line); ^/rj^ = 1, a = 2. ^3/2/^^ reaches a plateau at 
Tr ^ oc "1713/2/ y/fri^. The band denotes the WMAP result for cold dark matter with a 
2(7 error. 
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Figure 2: Relic gravitino density for different values of reheating temperature and grav- 
itino mass. ^Irf = 1. rrig = 1 TeV, which implies 777.3/2 < 0.1 TeV for a stable gravitino. 
For Tji> T^, 03/2/i^ is independent of Tr and m3/2- 
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One may also use eq. (35) to determine the range of gluino masses consistent with 
the WMAP result for cold dark matter. Varying QcBuh^ and /(q,)^(T*) in the ranges 
specified above we find, 

/ 2\l/6 

mp = (0.5...2.0) TeV (yj . (36) 

Hence, the hypothesis that gravitinos are the dominant component of dark matter will 
be tested at LHC ! 

The range for the gluino mass given in eq. (36) has been obtained in the case of 
gaugino mediation where 7713/2 = {2ri/gl) {M/Mp)mg (cf. (5) and (6)), with 77 = 0{1). 
For gravity mediation [10], one obtains the same results with rj replaced hy rj' — rj Mp/M. 
m3/2 and rrig now have the same order of magnitude, but the gravitino can be the LSP 
without fine tuning. The range for the gluino mass remains unchanged unless M is 
smaller than Mp by several orders of magnitude. In the case of gauge mediation [11], 
77 has to be replaced by rj' — r]S7r'^{X)/M where {X) is the messenger scale. The mass 
range (36) for the gluino mass is then obtained if M is of order Stt^ {X) . Note that 
the rapid decrease of gauge couplings at high temperature occurs independently of the 
supersymmetry breaking mechanism. 

Our results have important consequences for leptogenesis [12] where the typical 
baryogenesis temperature is Tb — O(10^° GeV) or larger [13]. According to previous 
studies this implies that unstable gravitinos have to be heavier than a few TeV [14, 15]. 
Stable gravitinos may have masses below 0{1 keV) [16] so that their mass density is be- 
low the critical density even when they are thermalized. Further, it has been shown that 
also gravitino masses 777,3/2 ~ 10 . . . 100 GeV can be consistent, which then constrains 
masses and couplings of other neutralinos and sleptons [17, 18] . 

Our analysis shows that there is no constraint on the reheating temperature for 
gluino masses below 0{1 TeV) (cf. fig. 2) if the gravitino is the hghtest supersymmetric 
particle. For rrig — 0{1 TeV) and reheating temperatures Tr > we find ^^3/2^1^ — 
QcDM^^ — 0.1, independently of 1713/2- 

The maximal value of the critical temperature is obtained for M ~ Mp and 
1713/2 ~ 100 GeV, so that the gravitino can still be the LSP for a gluino mass 0{1 TeV). 
This yields T™^ ~ 10^° GeV, which happens to coincide with the typical leptogenesis 
temperature. Hence, for a reheating temperature Tr larger than the leptogenesis temper- 
ature Tb, relic gravitinos always have the observed dark matter energy density Qcdm^^ 
if the gluino mass is 0(1 TeV). In this way the supersymmetry breaking scale in the ob- 
servable sector is directly determined by the dark matter density flcBMh^, independently 
of the supersymmetry breaking scale in the hidden sector! 

The interplay of particle physics and cosmology relates some properties of the uni- 
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verse to properties of elementary particles. Of particular interest is the composition of 
the present energy density flh^. Leptogenesis explains the baryon density flbh^ in terms 
of neutrino masses and mixings. As we have seen, for stable gravitinos the dark matter 
density flcBuh'^ is then determined by the gluino mass, i.e., the supersymmetry breaking 
scale in the observable sector, which may also be responsible for the dominant contribu- 
tion to the cosmological constant. 
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